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Abstract 

We study the potential of a very large volume underwater Mediterranean neu- 
trino telescope to observe neutrinos from supernova (SN) explosions within 
our galaxy. The intense neutrino burst emitted in a SN explosion results 
in a large number of MeV neutrinos inside the instrumented volume of the 
neutrino telescope that can be detected (mainly) via the reaction v e + p — > 
e + + n. In this study we simulated the response of the underwater neutrino 
telescope to the electron antineutrino flux predicted by the Garching model 
for SN explosions. We assumed that the neutrino telescope comprises 6160 
direction sensitive optical modules, each containing 31 small photomultiplier 
tubes. Multiple coincidences between the photomultiplier tubes of the same 
optical module are utilized to suppress the noise produced by 40 K radioactive 
decays and to establish a statistical significant signature of the SN explosion. 
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1. Introduction 

Core collapse supernovae produce huge bursts of MeV neutrinos which 
carry most of the gravitational energy (~ 99%) of a massive star in a 
timescale of a few seconds. The released energy (~ 10 53 ergs) is expected to 
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be roughly equipartioned among the different flavors of neutrino^] emerging 
promptly from the dense core, a few hours before the optical signal from 
the stellar envelope. The neutrino burst's time, flavor and energy struc- 
ture carry information about the explosion mechanism |l( an d they can also 
be used to study vacuum and matter neutrino oscillations the neutrino 
mass hierarchy [3[ or even exotic physics phenomena related to large extra 
dimensions and couplings to axions [4|. Since the SN1987A explosion when 
24 neutrinos were detected by the Kamiokande-II, 1MB and Baksan, many 
other experiments were designed [a, @| to observe and study SN neutrinos. 
These detectors can also provide an early alarm for a SN explosion offering 
the chance to the astronomical community to make unprecedented observa- 
tions of the early turn-on of the supernova light curve. For this purpose a 
SuperNova Early Warning System (SNEWS |7|) already exists providing also 
a trigger to detector facilities which are not able to trigger by themselves on 
a supernova signal. 

At the time of this writing, supernova neutrinos are the only neutrinos 
detected beyond our solar system. However, a lot of effort has been invested 
in observing high energy neutrinos from galactic and extragalactic sources. 
Very large volume neutrino telescopes are designed [sj and are already oper- 
ating ||. Such telescopes are consisting of thousands Photomultiplier Tubes 
(PMT) and instrumented volumes of the order of km 3 . Neutrino interaction 
products, mainly muons and electrons, are detected through the Cherenkov 
light produced in water or in ice, providing a neutrino detection energy- 



threshold of about 100 GeV. However, it has been shown [10] , I 111 . 1 121 . 113L 114 
that bursts of low energy (MeV) SN neutrinos are detectable in neutrino tele- 
scopes in ice, by observing a global rise of the PMTs counting rates which 
exceeds significantly the background rate. Since undersea neutrino detectors 
are suffering from background noise, due to the 40 K natural radioactivity of 
the sea-water, it was generally believed that only very dense detector con- 



figurations could observe SN signatures [15|. In this paper we demonstrate 



that deep undersea neutrino telescopes equipped with optical sensors that 
comprise many small PMTs, are capable of detecting burst of SN neutrinos 
without imposing any restrictions on the detector architecture. 



1 About 1% of this energy is emmited in the form of electron neutrinos from an initial 
neutrinization burst lasting a few ms. 
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2. The Underwater Mediterranean i/-Telescope 



In this study, we assume an underwater neutrino telescope following the 
Mediterranean KM3NeT architecture [gj. The KM3NeT detector will con- 
sist of several hundreds of vertical structures (Detection Units - DUs), which 
carry photo-sensors and devices for calibration and environmental measure- 
ments, arranged vertically on Storeys. Each Storey will support two photo- 
sensors. The photo-sensor unit is a digital optical module (DOM) consisting 
of a 17 inch diameter pressure resistant glass sphere, housing 31 3-inch small 
photomultiplier (small-PMT) tubes, their high- voltage bases and their inter- 



faces to the data acquisition system with nanosecond timing precision 16 
The segmentation of the photocathode area in such a Multi-PMT DOM aids 
in distinguishing single-photon from multi-photon hits, providing thus a bet- 
ter background rejection and trigger efficiency The front-end electronics will 
utilize time over threshold (TOT) measurements that carry amplitude and 
timing information, which allow for reconstruction of the original signal. The 
readout system will support the broadcasting of all the digitized data to the 
shore. 

In this work, we studied the SN discovery potential of such a //-Telescope 
configuration, deployed in a Mediterranean site of 3500 m depth, consisting 
of 154 DUs, with a mean horizontal separation of 180 m. Each DU was 
arranged in 20 storeys (6 m long bars), with a vertical distance between 
them of 40 m. Each storey carries two photo-sensors, one at each end of the 
bar. 



3. Discovery Potential for SN explosions 

Assuming that the expected number of background events in a time in- 
terval T is /Zbck(^) and that m(T) events have been observed in this time 
period, the condition that there is less than 5.72 x 10~ 7 probability the ob- 
served events to be produced solely by the fluctuating background is: 

i-£iWZ2f <i-«r/(4) (i) 

i=0 

In case that m(T) is large, the above formula can be replaced by its Gaussian 
approximation, s(T) > 5^/yUbck, where s(T) = m{T) — fihck(T) is the excess 
of events with respect to the expected background. Obviously, the aim of 
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this study is to define an analysis strategy in order to maximize the following 
quantity: 



(2) 



In the case that the excess s(T) is due to neutrinos from a SN explosion, 
it will depend on the proximity of the SN to earth, i.e. it will be inverse 
proportional to the squared distance (d 2 ). Consequently, the 5a Confidence 
Level (CL) requirement in ([1]) determines the maximum distance of a typical 
SN in order to be detected, or equivalently, given the distance of the SN, the 
statistical significance of the observation can be evaluated. 

4. SN model and neutrino interactions 

The SN explosion model used in this work, is the so called "Garching 



Model" [1?J which assumes an initial star with mass of 8-10 solar masses, 
which represent 30% of all the studied SNe. 

Among the interactions of these SN neutrinos in water, inside the in- 
strumented volume, we consider only the inverse beta decay (IBD) of the 
electron antineutrino {v e + p — > e + + n) which is the most dominant pro- 
cess^. We used the IBD total cross section a(E Pe ) as well as the differential 
cross sections 

aE e + dcosu 
(where 9 is the angular difference between the positron direction and the 



neutrino direction) from 18], which are very accurate for the energy range 
of SN neutrino^. 

The time-integrated electron antineutrino differential flux (integrated for 
an observation time interval, T, which starts with the bounce) is expressed 
as: 

;f{Ep e ;{E Pe },a)dt 



dE Pe And 2 J (E Pe )(ty 

where L Ve and E Ve are the luminosity and the energy of the antineutrino re- 
spectively, whilst d is the distance of the SN from the earth. f(Ep e ; (Ep e ), a) 



2 The neutrino energy threshold for the IBD process is 1.8 MeV, accounting for about 
95% of the signal. 

3 NNLO corrections are less than 1%. 
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is the energy spectrum of the electron antineutrino as parametrized in [19 
where (E Pe ) is the mean energy of the neutrino and a a parameter represent- 
ing the amount of spectral pinching. The time evolution of L Pe , (E Pe ) and a 



are given by the Garching Model. The integration time T is chosen as in [14 
in order to maximize the discovery potential for SN detection (i.e. Eq. (j2J)) 
and it is found to be T = 2.95 s, using the time structure of the neutrino 
burst predicted by the Garching Model. 

5. Simulation studies 

The positron emitted in SN neutrino induced IBD process is sufficiently 
energetic (mean energy 17.19 MeV) to produce Cherenkov radiation in the 
sea water, which is detectable by the photo-sensors of a neutrino telescope. 
The response of the detector to the IBD positrons and background processes 



is simulated by utilizing the HOURS [20[ package (Hellenic Open Univer- 
sity Reconstruction and Simulation package), which describes in detail all 
the relevant physical processes (e.g. secondary interactions, propagation of 
Cherenkov photons), the response of the PMTs and the digitization electron- 
ics. Using the simulated events we evaluated the effective volume, V e g, for 
IBD interactions and the interaction density, Pmt{d\T), which expresses the 
number of IBD neutrino interactions per unit volume, occurring inside an 
observation time interval T, caused by the explosion of a SN at distance d. 
The effective volume is estimated as the water volume, around a DOM, in 
which we have simulated IBD positron production, weighted by the fraction 
of the interactions observed by the DOM, under certain observation criteria. 
The interaction density is found by estimating the integral, 

f°° d$>„ 

p int {d;T) = p targct / —-^a{E Pe )dE 9e 

J Ve 

by Monte Carlo integration, where ptarget is the free proton density in sea 
water at the detector deployment depth and the cross sections are taken as 
described in Section 4. As an example, for a SN explosion at a distance of 
d — 10 kpc and for T = 2.95 s the interaction density at a depth of 3500 m is 
estimated to be p int = 0.078 m~ 3 . Finally the expected number of observed 
IBD interactions, from a SN explosion at a distance d is evaluated as, 

AU = Pmt (d;T) x Kff (3) 
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Figure 1: The expected number of observed events (active DOMs), in an observation time 
interval T — 2.95 s, as a function of the PMT multiplicity for three type of events: 40 if 
decays, atmospheric muons and SN neutrinos. The SN explosion is assumed to occur at 
a distance of 10 kpc and the i/- Telescope is assumed to be deployed at a depth of 3500 m 
containing 6160 DOMs. 
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As already mentioned, the main background source to SN neutrinos is 
the natural radioactivity of 40 K in the sea water. The 40 K decay rate is 
about 13.6 x 10 3 m~ 3 s -1 , five orders of magnitude greater than the expected 
IBD rate from a SN at 10 kpc. However, the released 40 K decay energy 
is relatively smalfl with respect to the energy of the SN neutrino induced 
IBD positron. The effective volume corresponding to the 40 K decays, is thus 
much smaller than the effective volume for SN neutrino IBD interactions. 
Furthermore the products of 40 K decays have less probability to activate 
synchronously several small PMTs of the same DOM. We define thus, as the 
main observation criterion, the PMT multiplicity, which is the number of 
PMTs in the same DOM that are simultaneously active within a coincidence 
time window of 100 ns. Fig. [T] presents the total number of active DOMs due 
to 40 K decays, collected during an observation time interval of T = 2.95 s, 
as a function of the PMT multiplicity, in comparison with the corresponding 
number due to SN neutrino induced IBDs, assuming that the SN explosion 
occurred at a distance of 10 kpc. Although the background due to 40 K 
decays is significantly reduced at high PMT multiplicities, it is expected that 
atmospheric muons will still contribute to the background up to very high 
PMT multiplicities, due to the causal connection between hits produced by 
muons. In order to quantify this muonic background we have used HOURS 



with generation input from the MuPage package [21| with an energy threshold 
as low as 50 GeV. The contribution of this type of background is also shown 
as a function of the PMT multiplicity in Fig. [TJ 

6. Observation Strategy and Results 

As shown in Fig. [TJ a high PMT multiplicity rejects the m K background. 
A further suppression of this background can be achieved when we apply 
selection criteria, based on the amplitude of the small-PMT signals. Fig. [2] 
(a,b and c) shows the probability of a DOM, activated by a 40 K decay or 
SN neutrino induced IBD, to exhibit a certain PMT multiplicity, for three 
different small-PMT pulse- amplitude thresholds, Qthres- As an example, in 
Fig. |2^a) a small PMT is considered active when its pulse amplitude ex- 
ceeds an amplitude-threshold corresponding to the average pulse amplitude 
produced by one photoelectron (pe), Qthres = 1 pe. The probability of the 



4 the mean electron energy is 1.27 MeV and the energy of the photon due to electron 
capture is about 1.46 MeV 
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Figure 2: The probability of a DOM to exhibit a certain PMT multiplicity for three 
different small-PMT pulse-amplitude thresholds: Qthres = 1 pe (a), Qthrcs = 2 pes (b) 
and Qthres = 3 pes (c). (d): The probability that a DOM activated by an atmospheric 
muon, with PMT multiplicities greater that 5, has simultaneously (within a time window 
of 100 ns) N active neighboring DOMs, with multiplicity greater than 2. 
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40 K background to produce high PMT multiplicity events decreases faster 
than the signal, as the threshold requirement is getting higher. However we 
found that this requirement has not any significant effect on the background 
produced by atmospheric muons. On the other hand, it is very improbable 
that the 40 K background and the SN neutrino induced IBDs will activate 
small PMTs in different DOMs (especially when high PMT multiplicity is 
required). On the contrary, muon tracks produce very often time-correlated 
hits in different DOMs, with high PMT multiplicity. This is demonstrated 
in Fig. Efd), where it is shown the probability that a DOM activated by an 
atmospheric muon, with PMT multiplicity greater that 5, has simultaneously 
(within a time window of 100 ns) N active neighbor DOMs, with PMT mul- 
tiplicity greater than 2. Rejecting thus all the active DOMs with PMT multi- 
plicities greater than 5, when neighboring with at least one, simultaneously- 
active DOM, with PMT multiplicity at least 3, the atmospheric muon back- 
ground is reduced by 94% whilst only 10% of the signal is lost. We used fully 
simulated signal (10 SN neutrino induced IBDs) and background (10 11 40 K 
and 10 6 atmospheric muon) events in order to estimate (using Eq. (jHJ)) the 
number of signal events from a SN explosion at 10 kpc, as well as the ex- 
pected number of 40 K and atmospheric muon background events, for various 
sets of event selection criteria. In the simulation we have assumed that the 
response function of the PMTs produced by a single pe is described by the 
function 

f(t) = at b e- L c 

that the charge resolution is 35%, the time jitter is 1 ns and the rise time of 
the pulse is 2.5 ns. We have also assumed that the digitization electronics 
provide the arrival time and the TOT of the pulse at the 30% of the mean 
pulse amplitude corresponding to one pe. In Fig. [3] it is shown the mean 
value of the PMT-waveform Time over Threshold duration as a function of 
the corresponding number of photoelectrons, for three type of events (IBD 
interactions, 40 K decays and atmospheric muon events). As an example, the 
condition TOT > 7.5 ns selects pulses corresponding to more than 1 pe. 

The set of the optimal selection criteria, which are maximizing the sensi- 
tivity of the detector as expressed in Eq. (j2J) are tabulated in Table [TJ With 
these event selection criteria we expect 118.7 signal events from a SN explo- 
sion at 10 kpc on top of 0.8 m K and 15.4 atmospheric muon background 
events. For this expected background, an observation with 5cr significance 
requires, according to Eq. (pQ), at least 39 total observed events (22.8 signal 
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Figure 3: The mean value of the TOT value with respect to the number of the photo- 
electrons for three type of events: IBD interactions, 40 K decays and atmospheric muon 
events. The error bars correspond to the RMS of the TOT distribution. 
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Table 1: Selection criteria that maximize the sensitivity of the detector. 



Description 


Optimum Value 


PMT multiplicity > 


5 


TOTthrcshoid > 


7.5 ns 


Number of PMTs (with TOT > TOT thrcs ) > 


2 


Neighbor DOMs (with PMT multiplicity >2) < 


1 



events). Consequently, the maximum distance of a SN explosion that could 
be observed by this detector is estimated, by scaling proportionally to the 
squared distance, to be (10 kpc x ^118.7/22.8) 22.8 kpc. 
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